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1.  Introduction. 

The  ability  of  energetic  charged  cosmic-ray  particles  to  arrive  at  a  par¬ 
ticular  location  within  the  magnetosphere  can  be  characterised  by  the  cutoff 
rigidity  for  the  location  and  the  direction  of  particle  arrival.  Knowledge  of 
the  cutoff  rigidity  is  a  pre-requisite  for  the  determination  of  the  total  radi¬ 
ation  dose  from  such  particles  at  the  location  from  the  direction  under  consider¬ 
ation.  Calculation  of  cutoffs  is  a  time  consuming  task,  which  can  be  expedited 
by  prior  knowledge  of  the  directions  from  which  primary  particles  of  any  rigidity 
are  able  to  reach  a  specific  location. 

The  general  aim  of  the  research  reported  here  was  to  extend  the  information 
previously  available  on  the  access  of  cosmic  rays  to  Earth  orbiting  spacecraft, 
at  altitudes  between  a  few  hundred  and  a  few  thousand  kilometres.  This  range 
encompasses  the  locations  encountered  by  all  but  the  lowest  orbit  satellites 
having  periods  of  a  few  hours  or  less.  There  were  several  specific  aims  within 
the  general  heading.  Particular  attention  was  to  be  paid  to  below-horizon 
directions  of  arrival;  to  estimates  of  the  effect  due  to  the  proximity  of  the 
solid  Earth;  and  to  attempts  to  identify  locations  and  altitudes  at  or  near 
which  interpolation  between  cutoffs  calculated  for  adjacent  locations  and/or 
directions  of  arrival  does  not  produce  usefully  accurate  results.  It  was  also 
planned  to  start  to  investigate  the  sensitivity  of  the  calculated  results  to  the 
various  available  models  of  the  geomagnetic  field,  and  to  assumptions  made  in 
connection  with  the  effect  of  the  atmosphere. 


2.  The  effect  of  altitude  on  access  from  below-horizon  directions. 

Directions  from  which  at  least  some  particles  are  able  to  arrive  at  a  given 
location  may  in  principle  be  determined  by  the  use  of  a  brute  force  method  simi¬ 
lar  to  that  needed  for  explicit  calculation  of  cutoffs.  A  very  large  comput¬ 
ational  effort  is  required  to  produce  the  necessary  5-dimensional  matrix  of 
values  using  the  standard  trajectory-tracing  technique.  The  variables  involved 
are  the  location  (latitude,  longitude,  and  altitude),  and  the  particle's  direct¬ 
ion  of  arrival,  expressed  in  terras  of  azimuth  and  zenith  angles.  It  was,  there¬ 
fore,  never  intended  that  the  project  should  be  undertaken  in  that  manner. 

An  alternative  way  to  approach  such  a  survey  is  to  consider  a  particular 
location  and  azimuth  of  arrival,  and  to  then  search  for  the  largest  possible 
zenith  angle  at  which  the  location  may  be  reached  by  an  incoming  primary  cosmic 
ray  particle.  A  technique  was  devised  to  search  rigidity/zenith-angle  space 
for  accessible  arrival  directions  (Humble  et  al,  1983,  attached  as  Appendix  B). 
The  search  proceeds,  for  a  given  location  and  azimuth,  in  the  directions  of 
increasing  zenith  angle  and  decreasing  rigidity,  until  no  further  accessible 
directions  of  arrival  can  be  found. 

Each  search  was  started  at  a  zenith  angle  of  100°,  above  the  local 
satellite/Earth  horizon  at  all  altitudes  considered  (Humble,  1983,  attached  as 
Appendix  B;  Table  1).  It  was  found  that,  for  zenith  angles  larger  than  that  of 
the  horizon,  primary  particles  are  only  able  to  reach  the  satellite  from  gener¬ 
ally  westerly  directions.  Such  particles  experience  a  V_  x  force  having  a  pos¬ 
itive  radial  component  in  the  final  stages  of  their  approach  to  the  satellite. 
Their  trajectories  consequently  have  positive  upwards  curvature,  and  the  local 
zenith  angle  of  arrival  can  be  larger  than  that  of  the  horizon.  The  range  of 
accessible  azimuths  increases  with  altitude,  as  would  be  expected. 
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In  general,  It  would  be  expected  that  the  largest  accessible  zenith 
angle  for  a  given  latitude,  longitude,  and  azimuth,  would  also  increase  with 
increasing  altitude.  This  was  found  to  be  the  case.  Defining  Z  to  denote  the 
largest  accessible  zenith  angle  found  at  any  altitude,  regardless  of  latitude, 
longitude,  or  azimuth,  the  following  results  were  obtained: 


Altitude 

Z 

400  km 

150 

600 

156 

800 

166 

1000 

172 

1250 

180 

More  detailed  results  are  given  in  the  two  publications  which  have  so  far  arisen 
from  this  work  (Humble  et  al,  1983:  Appendix  A;  Humble,  1983:  Appendix  B). 

Note  that  a  Z  of  180°  for  a  particular  altitude  means  that  primary  particles 
are  able  to  reach  a  satellite  at  that  altitude  from  the  nadir  direction,  directly 
underneath  it,  at  at  least  one  location  somewhere  on  the  orbit.  This  does  not, 
however,  mean  that  all  possible  directions  of  arrival  are  accessible  at  such  a 
location.  Easterly  arrival  directions  at  below-horizon  zeniths  are  still 
forbidden. 

As  the  altitude  increases  above  1250  km  a  range  of  large  zenith  angles  in 
generally  easterly  directions  begins  to  become  accessible,  whilst  smaller  (but 
still  below-horizon)  zenith  angles  in  the  same  azimuthal  directions  remain  in¬ 
accessible.  The  higher  altitude  permits  particles  approaching  the  general  region 
of  the  satellite  at  low  altitudes  from  the  west  to  pass  above  the  top  of  the 
atmosphere  more  than  one  gyro-radius  beneath  the  satellite.  Such  particles  will 
arrive  at  the  satellite  from  easterly  directions  at  large  zenith  angles. 

The  phrase  "generally  easterly"  used  above  requires  qualification.  The 
precise  azimuths  at  which  these  effects  occur  are  latitude  dependent.  The 
largest  zenith  angles  of  access  occur  in  directions  which  are  slightly  equator- 
ward  of  west  in  both  hemispheres.  The  curve-under  effect  occurs  in  the  set  of 
opposite  directions,  slightly  to  the  poleward  of  east  in  both  hemispheres. 

Preliminary  surveys  were  conducted  at  altitudes  of  2400  and  3600  km.  These 
computations,  like  all  others  in  this  study,  were  performed  for  20°  intervals  in 
latitude  and  60°  intervals  in  longitude.  At  2400  km  the  curve-under  phenomenum 
was  found  to  exist  only  at  latitudes  between  20°N  and  20°S.  At  these  latitudes 
the  effect  is  also  markedly  longitude  dependent,  due  to  the  asymmetric  nature  of 
the  geomagnetic  field.  Accessible  easterly  zenith  angles  were  found  to  commence 
at  the  nadir,  zenith  angle  180Q,  and  to  extend  to  zenith  angles  of  about  140°, 
with  a  couple  of  cases  running  up  to  110°,  well  above  the  satellite/Earth  horizon 
at  this  altitude. 

Similar  results  were  obtained  for  3600  km  altitude,  except  that  the  range 
of  allowed  easterly  zeniths  is  much  larger  and  the  latitude  range  extends  to  40° 
in  both  hemispheres.  The  longitude  effect  is  not  so  noticeable.  Furthermore, 
a  considerable  number  of  cases  were  found  at  this  altitude  for  which  some  east¬ 
erly  azimuths  are  accessible  at  all  zenith  angles.  In  these  cases  the  range  of 
zeniths  accessible  by  particles  from  the  west  looping  below  the  satellite  over¬ 
laps  the  range  accessible  by  particles  of  other  rigidities  approaching  more 
directly  from  the  east.  Due  to  the  nature  of  the  search  algorithm  employed, 
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it  is  suspected  that  a  few  more  such  cases  may  exist  than  have  actually  been 
detected.  Detailed  checking  in  a  few  directions  suggests  that  the  algorithm 
probably  detects  at  least  90%  of  the  occurrences. 


3.  The  Interpolation  Problem. 

The  results  obtained  in  the  course  of  the  above  calculations  show  that  there 
are  some  regions  and  directions  at  which  simple  interpolation  of  cutoffs  between 
adjacent  directions  and/or  locations  is  able  to  produce  usefully  accurate  results, 
whilst  at  other  directions/locations  this  is  not  so.  In  general,  cutoff  inter¬ 
polation  in  below-horizon  directions  is  practicable  at  those  directions  and 
locations  at  which  the  gradient  of  the  maximum  accessible  zenith  vs.  azimuth 
curve  (Humble, 1983:  Appendix  B;  Figure  1)  is  moderate.  Useful  interpolation  is 
impracticable  wherever  that  curve  is  steep,  since  discontinuities  in  maximum 
accessible  zenith  angle  occur  quite  frequently  in  that  situation.  The  dis¬ 
continuities  are  due  to  the  increasing  or  decreasing  of  the  influence  of 
particular  components  of  the  geomagnetic  field  as  azimuth  of  arrival  changes. 

The  accuracy  of  the  interpolation  process  varies  both  with  altitude  and 
with  azimuthal  direction.  The  most  reliable  values  can  be  obtained  between 
azimuths  of  195°  and  345°  at  3600  km  altitude,  reducing  to  a  rather  smaller 
azimuthal  range,  of  about  90°,  at  400  km.  This  range  is  centred  about  dir¬ 
ections  which  are  approximately  (270  -  latitude)0,  significantly  to  the 
equatorward  of  west  at  all  but  the  lowest  latitudes. 


4.  Comparison  of  Geomagnetic  Field  Models. 

The  present  project  is  based  on,  and  has  made  use  of  some  results  from, 
work  performed  previously.  That  work  used  the  then  best  available  numerical 
model  of  the  internal  geomagnetic  field,  being  an  8th  order  model  based  on  a 
1975  field  model  extrapolated  forward  to  epoch  1980.0  by  use  of  associated 
secular  drift  coefficients.  It  was  known  at  that  time  that  this  model  was  not 
ideal,  and  that  more  accurate  models  would  shortly  become  available. 


About  the  time  the  present  project  commenced  a  10th  order  definitive 
International  Geomagnetic  Reference  Field  for  epoch  1980.0  was  published. 
Accessibility  calculations  were  immediately  performed  for  selected  locations 
using  this  field  model,  in  order  to  assess  the  sensitivity  of  the  results  to 
the  choice  of  field,  and  to  determine  if  it  would  be  necessary  to  recalculate 
any  or  all  of  the  results  previously  obtained.  The  10th  order  field  model 
requires  the  use  of  approximately  25%  more  computer  time  than  do  the  8th 
order  models,  and  comparative  calculations  were  therefore  also  rierformed 
using  the  10th  order  field  arbitrarily  truncated  to  8th  order.  The  sensitivity 
of  the  results  to  the  choice  of  field  model  was  found  to  be  location  dependent, 
being  most  noticeable  in  the  South  Atlantic  and  South  American  regions  f  & 

(Humble  et  al,  1984).  (This  paper,  which  is  attached  as  appendix  C,  reports 
an  amalgamation  of  work  carried  out  under  the  present  contract  and  related  worUi/ 
performed  separately  by  AFGL  staff).  The  sensitivity  in  other  regions  is  not  so 


5.  The  effect  of  the  atmosphere. 

This  aspect  of  the  project  was  directed  towards  determining  the  effect  on 
the  results  of  the  use  of  a  realistic  model  atmosphere  to  determine  cosmic-ray 
absorption,  rather  than  making  use  of  the  standard  assumption  that  the  atmos¬ 
phere  is  completely  opaque  to  any  primary  particle  which  descends  below  an 
altitude  of  30  km.  Only  very  preliminary  calculations  were  completed,  amounting 
to  little  more  than  the  necessary  checking  of  the  computer  programmes  involved. 
The  few  available  results  appear  to  indicate  that  the  30  km  limit  is  a  sensible 
approximation.  This  result  had  been  anticipated;  however,  considerably  more 
work  will  be  required  before  any  definitive  conclusions  can  be  drawn  from  this 
aspect  of  the  project. 
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COSMIC  RAY  ACCRSS  TO  .SATELLITES  FROM  LARGE  ZENITH  ANGLES 

J.E. Humble 

Department  of  Physics 
University  of  Tasmania 
Hobart,  Tasmania 
Australia  7001 
and 

D.r. Smart  and  M. A. Shea 
Air  Force  Geophysics  Laboratory  (PUG) 

Hanscom  AFR ,  MA  01731 
U.S.A. 

ABSTRACT 

Ac  tlic  altitude  of  cartli-orbiting  satellites  increases 
it  becomes  progressively  easier  for  primary  cosmic- 
ray  particles  to  gain  access  to  them  from  directions 
below  the  geometric  horizon.  Results  are  presented 
from  a  comprehensive  survey  of  the  largest  accessible 
zenith  angles  for  a  range  of  altitudes  and  geographic 
locations.  The  search  has  disclosed  that  primary 
particles  are  able  to  reach  a  satellite  at  an  altitude 
of  1250  km  from  zenith  angles  as  large  as  178°. 

1.  Introduction 

Some  interest  has  been  expressed  in  recent  times  in  the  ability  of 
primary  cosmic-ray  particles  to  roach  satellites  in  low  altitude  earth- 
orbit  from  directions  below  the  local  horizon.  At  the  Paris  Cosmic  Ray 
Conference,  Humble  et  al.  (1981)  showed  that  primary  particles  are  able 
to  reach  a  satellite  at  a  zenith  angle  of  120°  from  a  range  of  western 
azimuths  at  all  the  geographic  locations  which  were  investigated.  The 
cal  dilations  were  substantiated  by  the  results  of  experiments  on  board 
the  11F.A0-C  satellite  (N.Lund,  private  communication  1981).  The 
investigation  reported  here  extends  these  studies  by  determining  in  a 
systematic  fashion  the  .largest  zenith  angles  accessible  to  primary 
particles  at  a  selected  range  of  azimuths  and  set  of  satellite 
altitudes  and  locations. 

Calculations  arc  reported  for  altitudes  of  400,  800  and  1250  km, 
for  latitudes  at  20°  intervals  from  40°N  to  40°S  at  each  altitude,  and 
longitudes  at  60°  intervals  starting  at  0°  at  each  latitude.  For  each 
location  thus  defined  a  suitable  range  of  azimuths  likely  to  be  accessible 
to  primary  particles  at  zenith  angles  of  120°  or  larger  (Humble  et  al., 
1981)  was  investigated. 

2.  Method 

The  calculations  have  been  carried  out  using  tbe  standard  traject¬ 
ory  calculation  program  (Shea  et  al.,1965).  For  continuity  with  earlier 
work  we  have  continued  using  the  former  International  Geomagnetic  Refer¬ 
ence  Field  for  1980.0,  as  extrapolated  from  ICRF  1975.0  with  the  use  of 
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secular  collection  coefficients  (JANA  Study  Group , 1 076) .  We  note  that 
future  calculations  should  use  the  new  definitive  ICKK  l'.'OO.O  model  result¬ 
ing.  from  the  Tdinburgh  meeting  of  1AGA  (l’eddie  ,  I  98?  )  . 

As  zenith  angles-  increase  at  a  particular  location  and  azimuth  of 
arrival  the  Horizon  Limited  Rigidity*  necessarily  decreases,  due  to  the  in¬ 
creased  curvature  required  for  the  final  stage  of  the  trajectory  (Humble  et 
a  1 .  , !  9  8 1  )  .  A  similar  situation  holds  true  for  particular  penumbra  1  bands. 

The  rather  simple  minded,  but  effective,  search  technique  used  was,  there¬ 
fore,  to  commence  calculations  at  zenith  angle  12(1°  and  a  large  rigidity, 
usually  around  30  GV.  If  this  trajectory  was  allowed  the  zenith  angle  was 
incremented  by  2°  whilst  if  it  was  forbidden  the  rigidity  was  decremented 
by  I?..  The  process  wan  then  repeated  at  the  new  zenith  or  rigidity. 

The  technique  is  fast,  in  that  the  majority  of  the  trajectories 
involved  are  simple  ones.  It  has  the  effect,  of  following  particular  pen¬ 
umbra  1  features  through  increasing  zenith  angles.  If  one  such  hand  ends 
as  zenith  angle  increases  the  next  lowest  (in  rigidity  space)  band  is 
automatically  found.  This  process  continues  until  the  Stonier  cone  is 
encountered.  The  technique  does  ignore  the  possibility  that  a  penumbra 1 
allowed  band  may  exist  at  relatively  high  rigidities  extending  to  larger 
zenith  angles  than  do  tho  lower  rigidity  bands  first  explored.  One  or  two 
such  examples  have  been  found,  the  clear  implication  being  that  great  care 
must  be  used  in  selecting  the  starting  rigidity. 

3.  Results  and  Discussion 

The  largest  accessible  zenith  angles,  Z,  found  at  each  location  are 
listed  below. 
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‘•,;Thc  highest  rigidity  with  which  primary  particles  can  reach  the  site. 
This  was  called  "Allowed  Rigidity"  in  our  1981  paper;  we  believe  the 
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Figure  1  shows  scans  of  the  largest  accessible  zenith  angles  found 
versus  azimuth,  for  one  site  at  each  latitude  at  400  km  altitude.  Note 
that  the  largest  zenith  angle  accessible  at  each  site  occurs  at  azimuths 
to  the  eejuatorward  of  West.  Tins  .is  true  in  nil  cases  which  we  have 
investigated.  Note  also  the  fine  details  plotted  for  20°N  300°.  Exten¬ 
sive  calculations  were  undertaken  at  this  site. 

The  distribution  of  Z  with  latitude  varies  according  to  the 
altitude  involved.  At  400  km  there  is  no  detectable  relationship  between 
Z  and  latitude.  However,  a  relationship  begins  to  be  apparent  at  800  km, 
and  is  much  more  obvious  at  1750  km,  at  which  altitude  the  largest  values 
of  Z  are  found  between  20°N  and  20°S.  The  Z  values  at  higher  lati¬ 
tudes  are  noticeably  reduced. 

The  largest  values  of  Z  found  at  each  altitude,  and  indicated  in 
the  table,  arc  148°  at  400  km,  166°  at  800  km,  and  178°  at  1250  km. 

These  compare  with  horizon  zenith  angles  of  110,  117,  and  123°  respective¬ 
ly.  The  listings  only  indicate  values  actually  found  at  a  selected  grid 
of  locations.  They  do  not  exclude  the  possibility  of  larger  values  being 
found  at  non-grid  locations.  The  limited  spot  checks  which  we  have  made 
at  400  km  only,  for  various  latitudes  between  40°N  and  20°S  have  revealed 
one  such  case,  an  allowed  trajectory  at  a  zenith  angle  of  150°.  Due  to 
the  method  used  wc  cannot  completely  guarantee  that  larger  angles  than 
those  listed  also  exist  at  the  grid  locations.  However,  we  believe  that 
significant  variations  are  unlikely. 

Tt  should  also  be  noted  that  the  calculations  are  artificial  in  one 
important  respect.  They  assume,  along  with  earlier  calculations,  that  al.l. 
particles  descending  below  a  local  altitude  of  30  km  undergo  interaction 
with  the  atmosphere.  The  general  agreement  referred  to  earlier  between 
the  HEA0-C  observations  and  early  calculations  for  400  km  altitude  suggest 
that  the  assumption  may  be  reasonable,  but  the  point  requires  elucidation. 

Acknowledgements 

The  research  at  the  University  of  Tasmania  was  supported  by  the  Air 
Force  Geophysics  Laboratory  by  grant  numbers  AFOSR-80-0232  and  AFOSR-82- 
0313. 

References 

Humble  J  E,  Smart  D  F  and  Shea  M  A  -  Cosmic  ray  access  to  earth  satellites 
from  bclow-horizon  directions.  17th  International  Cosmic  Ray 
Conference,  Conference  Papers  JjO,  270-273  (1981) 

IACA  Division  1 ,  Study  Group,  International  Geomagnetic  Reference  Field 
1975,  J.Geophys.Res.  ,  8J^,  5163  (1976) 

Peddle  N  W  -  The  International  Geomagnetic  Reference  Field:  the  third 
generation.  J. Geomagnetism  and  Gcoclectricity  (1982)  in  press 
Shea  M  A,  Smart  I)  F  and  McCracken  K  G  -  A  study  of  vertically  incident 

cosmic  ray  trajectories  in  a  sixth  degree  simulation  of  the  geomag¬ 
netic  field.  J.Geophys.Res.  70,  4117  (1965) 


Proc.  ASA  5  (2)  1983 


11 


Contributions  265 


A  —  (A//(o)underj/oun<l  —  0. 14  (A///<,)|>jeulron 

Where  A/  is  (he  deviation  of  the  daily  average  intensity  /  from 
a  reference  intensity  /..  The  latter  was  chosen  as  the  value  of 
daily  average  intensity  on  1 1  July,  one  day  before  the 
commencement  of  the  decrease. 

The  A  variations  are  shown  in  Figure  2(c).  rhey  exhibit  a 
prominent  27-day  periodicity,  clearly  not  initiated  by  the 
Forbush  Decrease.  It  seemed  that  the  variations  might  be  due 
to  the  interplanetary  North-South  anisotropy  for  which  firm 
evidence  has  emerged  during  the  last  decade  (Swinson,  1971). 
this  is  a  B  x  ?(/i)  phenomenon,  where  B  is  the  interplanetary 
magnetic  field  (IMF)  vector  and  V(n )  is  the  density  gradient 
vector  of  cosmic  ray  particles,  directed  radially  outward  from 
the  Sun.  The  resulting  particle  drift  normal  to  the  plane  of  the 
ecliptic  produces  a  diurnal  variation  of  intensity  (in  sidereal 
time)  and  a  daily  average  intensity  term.  It  is  the  latter,  known 
as  the  North-South  (N-S)  asymmetry,  that  is  of  interest  here. 

The  N-S  intensity  1(6),  at  a  latitude  6,  varies  as  sin  6.  Thus  it 
is  asymmetric  with  respect  to  the  hemisphere  of  observation, 
an  increase  in  one  hemisphere  being  mirrored  b>  '  :  crease  in 

the  other  and  vice  versa.  Moreover,  the  asymmetry  reverses 
direction  with  reversals  of  direction  of  the  sectorised  field 
vector  B.  It  has  a  flat,  -/%,  rigidity  spectrum.  The  magnitude 
of  the  asymmetry,  at  mid-latitudes  of  observation,  is  usually 
small,  being  -0.1%.  Thus  it  can  easily  be  swamped  by  large 
isotropic  variations. 

Until  now,  only  the  large  multi-directional  telescope  system 
at  Nagoya,  Japan,  has  been  able  to  detect  the  N-S  asymmetry 
with  sufficient  accuracy  for  day-to-day  variations  to  be 
observed.  Isotropic  effects  are  removed  by  recording  the 
difference  of  intensity  between  north-pointing  and  south- 
pointing  (equatorial  scan)  telescope  components.  The  daily 
difference  measurements,  known  as  GG,  have  been  observed  to 
follow  a  27-day  variation,  within  which  the  peaks  tend  to 
correspond  with  the  Towards  direction  of  the  sectorised  IMF 
and  the  troughs  with  the  Away  direction.  This  expected 
correspondence  has  been  achieved  on  more  than  TQ*h  of  daily 
occasions  over  several  years  of  observation  (Mori  and 
Nagashima  1979). 

In  Figure  2(c)  the  variations  of  GG,  inverted  because  of  the 
latitude  asymmetry,  are  compared  with  the  variationsof  A.  It 
is  clear  from  the  close  agreement  between  them  (linear 
correlation  coefficient  r-0.62)  that  the  systematic  variations 
of  GG  and  A  are  due  to  the  same  phenomenon.  That  they  are 
due  to  the  N-S  asymmetry  can  be  seen  from  the  striking 
correspondence  of  the  peaks  and  troughs  of  the  variations  with 
the  Away  and  Towards  polarities  of  the  IMF  (Figure  2(d) ). 

Conclusion 

The  most  promising  feature  of  these  initial  results,  achieved  by 
quite  different  methods  in  the  two  hemispheres,  is  the  fidelity 
with  which  the  variations  from  Nagoya  and  Mawson  have 
tracked  each  other,  not  only  over  the  period  shown,  but  in 
succeeding  months.  It  means  that  a  sufficiently  high  degree  of 
stability  of  observation  has  been  attained  as  to  enable,  for  the 
first  time,  various  kinds  of  long-term  trends  in  the  average 
intensity  underground  to  be  investigated  in  both  hemispheres 


and.  further,  differences  in  trends  between  the  hemispheres, 
for  which  there  is  already  some  evidence  in  the  observations 
reported  here. 
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On  Directions  from  which 
Cosmic  Rays  may  Reach 
Earth  Satellites 

J.  E.  Humble,  Physics  Department,  University  of  Tasmania 


Introduction 

it  is  easy  to  feel,  on  an  intuitive  basis,  that  primary  cosmic 
rays  of  at  least  some  energy  range  will  be  able  to  reach  a 
satellite  in  earth-orbit  from  all  directions  above  the  local 
geometric  horizon  and  from  essentially  no  directions  below  the 
horizon.  This  argument  suggests  that  directional  cosmic  ray 
detectors  on  board  satellites  may  safely  assume  that  the  only 
particles  received  from  earthward  directions  must  be  splash  or 
albedo  particles  from  the  atmosphere,  and  that  any  equipment 
liable  to  possible  radiation  damage  should  be  preferentially 
located  on  the  earthward  side  of  the  satellite.  In  actual  fact 
primary  cosmic  rays  can  have  quite  sharply  curved  trajectories 
in  the  magnetosphere,  thereby  enabling  some  particles  to  gain 
access  to  a  satellite  from  a  range  of  directions  below  the  local 
geometric  horizon.  The  present  work  presents  recent 
quantitative  results  on  this  matter. 

The  Calculations 

Calculations  have  been  carried  out  using  the  standard 
computer  program  for  computing  cosmic  ray  trajectories  in  the 
magnetosphere  (Shea  ef  of.  1965).  For  continuity  with  earlier 
work  the  former  International  Geomagnetic  Reference  Field 
for  1980.0,  as  extrapolated  from  1GRF  1975.0  by  use  of 
secular  correction  coefficients  (IAGA  Study  Group  1976)  has 
been  used.  The  differences  caused  by  the  use  of  this  model  of 
the  geomagnetic  field  rather  than  the  more  recently  published 
definitive  IGRF  1980.0  model  (Peddie  1982)  will  be  the  subject 
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of  further  investigation.  Preliminary  computations  suggest  that 
the  differences  will  not  be  large. 

The  procedure  adopted  is  to  calculate,  by  numerical 
integration  of  the  equations  of  motion,  the  trajectory  of  a 
negatively  charged  particle  leaving  the  satellite  in  a  specified 
direction.  This  trajectory  will  be  precisely  the  same  as,  but 
traversed  in  the  opposite  direction  to,  that  of  a  positively 
charged  particle  arriving  at  the  satellite  from  the  specified 
direction.  If  the  trajectory  of  the  negative  particle  enters  the 
atmosphere,  at  an  altitude  of  30  km  or  less  above  mean  sea- 
level,  it  is  unlikely  that  an  inbound  positive  particle  would 
emerge  from  the  atmosphere  in  that  particular  direction,  since 
it  would  have  been  destroyed  in  interactions  with  atmosphere 
particles  before  doing  so.  The  particular  direction  of  arrival  is, 
therefore,  considered  to  be  forbidden  to  incoming  particles. 
The  argument  is  stronger  yet  in  the  event  that  the  trajectory 
intersects  the  solid  earth. 

The  only  significant  force  acting  on  the  cosmic  ray  particles 
in  the  magnetosphere  is  v  x  B,  and  it  is  obvious  that  the  only 
particles  able  to  reach  the  satellite  from  directions  below  the 
local  horizon  are  those  for  which  v  x  B  has  a  positive  radial 
component  in  regions  close  to  the  satellite.  Broadly,  this  will 
tend  to  facilitate  access  from  below  horizon  zenith  angles  for 
particles  arriving  from  westerly  azimuths  (Humble  el  al.  1981), 
whilst  inhibiting  such  access  from  generally  easterly  directions. 

Results 

Calculations  of  this  type  produce  a  plethora  of  results 
encumbered  with  a  number  of  parameters.  The  specific 
variables  are  the  location  (latitude,  longitude  and  altitude)  of 
the  satellite  and  the  magnetic  rigidity  and  direction  of  arrival 
(expressed  in  terms  of  zenith  and  azimuth  angles)  of  the 
particle  whose  trajectory  is  being  investigated.  The  present 
calculations  have  been  directed  towards  finding  the  largest 
zenith  angles  with  which  primary  cosmic  ray  particles  can 
reach  the  satellite  at  each  of  a  range  of  locations.  Searches  for 
such  large  accessible  zenith  angles  have  been  carried  out  at 
locations  specified  by  all  possible  combinations  of  latitude  at 
20°  intervals  from  40°N  to  40°S  inclusive,  longitude  at  60° 
intervals  starting  at  0°,  and  altitudes  of  400,  800  and  1230  km 
above  the  mean  spherical  earth.  A  smaller  range  of  latitudes 
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The  variation  of  largest  accessible  zenith  angle  with  geographic 
location  at  400  ka  altitude. 
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has  also  been  investigated  at  altitudes  of  600  and  1000  km. 
Some  additional  azimuths  were  also  investigated  at  400  km 
altitude,  with  the  aim  of  detecting  any  fine  scale  dependence  of 
largest  accessible  zenith  angle  with  azimuth. 

At  each  location  searches  have  been  carried  out  at  30° 
azimuth  intervals,  ranging  larger  or  smaller  as  required  from 
233°.  These  generally  westerly  directions  are  those  previously 
mentioned  for  which  v  x  B  is  known  to  have  a  positive  radial 
component  in  the  region  of  the  satellite.  The  risk  that  the  use 
of  azimuths  so  pre-determined  might  lead  to  a  self-fulfilling 
prophecy  in  the  final  results  was  avoided  by  taking  care  to 
extend  most  of  the  searches  in  azimuth  until  it  was  evident  that 
no  further  possibility  of  access  at  large  zenith  angles  remained. 

The  largest  accessible  zenith  angles  found  at  each  altitude, 
regardless  of  satellite  geographic  position  or  of  particle 
azimuth  of  arrival,  are  shown  in  Table  1 .  It  can  be  seen  that 
the  increase  in  the  angle  with  altitude  is  essentially  linear,  to 
within  the  coarseness  (2°)  of  the  zenith  steps  used. 

The  programme  used  to  obtain  the  above  results  makes  a 
complete  set  of  calculations  for  any  specified  location  and 
azimuth  of  particle  arrival.  For  each  set  of  initial  conditions  it 
searches  rigidity-zenith  space  for  the  largest  zenith  angle  for 
which  primary  cosmic  ray  access  can  be  confirmed.  The 
algorithm  used  (Humble  el  al.  1983)  is  satisfied  by  finding  a 
single  allowed  trajectory  at  a  further,  larger,  zenith. 
Consequently,  the  rigidity  of  the  particle  for  which  an  allowed 
trajectory  is  found  is  only  a  pseudo  randomly  chosen  member 
of  the  set  of  rigidities  with  which  particles  may  be  able  to 
arrive  at  the  satellite  from  the  direction  under  consideration.  It 
is  necessary  to  recognise  this  restriction  when  considering  the 
rigidities  quoted  in  the  final  column  of  Table  I . 

There  are  considerable  variations  between  the  largest 
accessible  zenith  angles  found  for  individual  longitudes  at  any 
given  latitude  and  altitude.  The  results  in  Table  2  show  this 
situation  for  an  altitude  of  400  km. 

The  inference  from  Table  2,  and  from  similar  computations 
at  other  altitudes,  is  that  the  largest  accessible  zenith  angle  is 
not  a  strong  function  of  latitude,  at  least  between  40°N  and 
40°S. 

The  azimuths  in  which  the  largest  accessible  zenith  angles 
occur  are,  however,  a  strong  function  of  latitude,  with  a 
marked  tendency  for  the  largest  zenith  angles  to  be  found  at 
azimuths  to  the  equatorward  of  west,  as  shown  in  Figure  I. 
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Figure  I.  The  variation  with  latitude  and  azimuth  of  largest  accessible 
zenith  angle  for  a  satellite  orbiting  at  400  km  altitude.  The  longitudes 
chosen  are  those  for  which  the  accessible  zenith  angle  was  largest  for 
each  latitude. 


Conclusion 

Satellites  in  earth  orbit  at  medium  altitudes  with  orbital 
inclination  of  40°  or  less  may  easily  be  reached  by  primary 
cosmic  ray  particles  from  zenith  angles  well  below  that  of  the 
local  geometric  horizon.  At  altitudes  of  12S0  km  and  above 
particles  can  reach  the  satellite  from  all  possible  zenith  angles. 
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Introduction 

After  21  years  the  Parkes  radio  telescope  has  undergone  a 
major  refit.  A  new  VAX-1 1/750  running  VMS  has  replaced 
the  aged  Digital  Equipment  Corporation  PDP-9  computer. 
Other  new  equipment  includes  two  systems  not  previously 
available  —  the  Mk  II  RING  communication  system  (Willing 
and  Abies  1983)  and  the  'observer  workstations’. 

A  new  observing  system,  called  HYDRA,  has  been 
developed  to  enable  the  general  user  to  make  the  best  use  of 
the  new  facilities. 

Specification  of  HYDRA 

Two  extreme  approaches  could  be  taken  when  designing  such  a 
system: 

(a)  The  piecemeal  approach  —  one  simple  program  working 
in  complete  isolation  for  each  type  of  observation. 

(b) .  The  ‘integrated’  approach,  in  which  the  hardware  and 

each  of  the  simple  programs  are  linked  into  a  complete 
environment  for  performing  astronomical  observations. 

It  is  neither  desirable  nor  practical  to  design  a  single 
monolithic  observing  program  to  control  the  whole  range  of 
observations  that  are  performed  at  Parkes.  A  continuum 
scanning  program  and  a  program  for  collecting  pulsar  data, 
for  example,  have  quite  different  requirements. 

The  design  of  HYDRA  allows  individual  components  to  be 
enhanced  or  replaced  with  ease.  To  this  end  we  have  provided 
a  software  framework  which  permits,  in  a  convenient  and 
standardized  manner,  communication  between  the  user,  his 
equipment,  his  data  files  and  the  telescope  control  and 
receiving  equipment.  This  framework  supports  software 
components  which  perform  the  data  collection  and  collation, 
as  well  as  post-observational  processing.  It  is  easy  to  add 
observing  programs  and  post-observational  processing 
programs  to  this  well-defined  framework  (as  long  as  a  number 
of  ground  rules  are  observed). 

The  following  requirements  were  taken  into  consideration 
during  the  design. 

(a)  Observers  often  wish  to  analyse  data  immediately;  they 
should  be  able  to  continue  observing  while  this  is  being 
done. 

(b)  Each  of  the  observing  programs  should  have  the  same 
user  interface.  This  interface  should  be  provided  via  a 
menu  system  in  which  all  of  the  data  input  to  the 
observing  program  are  specified.  The  menu  system  should 
be  many-layered,  allowing  the  user  to  control  the  various 
parts  of  the  system  separately.  A  help  facility  should  be 
provided  to  explain  the  function  of  that  layer,  and  to 
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ABSTRACT 


Cosmic  cay  vor!  ic:al  cutoff  rigidities  at  sea  Jev.-l  have  been  calculat¬ 
ed,  using  the  tra j ect ory -trac i ng  method,  for  a  number  of  different  epochs. 
These  calculation:;  have  been  carried  out  for  a  world-wide  grid  of  locations, 
and,  in  on  effort  to  locate  the  cosmic  ray  equator,  for  a  fine  grid  in  the 
equatorial  region.  Comparison  of  the  vertical  cutoff  rigidity  values  obtained 
using  the  international  Geomagnetic  Reference  Field  model  for  1930.0  with  those 


obtained  iron  previous  models  shows  systematic  significant  chang 


ion  iri  the 


Atlantic  Ocean  region  and  over  South  America.  The  differences  are  greater 
t'nor.  those  predicted  utilizing  the  older  field  models  with  their  predicted 
secular  change.  The  cutoff  rigidity  values  calculated  using  the  new  ISRF 
19iiJ.O  field  model  appear  to  be  in  better  agreement  with  data  from  cosmic  ray 
latitude  surveys  in  the  Atlantic  Ocean  region.  The  changes  in  the  cosmic  n.y 
equator  are  asymmetrical  with  essentially  no  changes  in  the  equator  position  in 
the  A&i  an  cur.i  Pacific  region,  but  with  significant  changes  in  the  South 
Anvr  i  can .  Atlantic  Ocean  and  West  Anican  regions.  Calculation:;  have  elro 
b-.-en  underi  akin  for  different  directions  of  arrival  for  a  satellite  ulitun  of 
<;urj  ;::n  aitihu.h-  using  the  predicted  i%0  field  model  and  1  he  inter,  in  1980  field, 
n'cdf'l  adopted  in  19 31 .  Some  ctif £‘  r  cnos  have  boon  found. 


1, _ TMVROiiHrTION 

The  cutoff  r Kii< lily  in  a  conrept  which  describes  the  magnet  j c  shieldii:  f 
provided  by  the  geomagnetic  field  against  the  arrival  of  charged  cosmic:  » 
par  tide- s.  Rigidity  in  a  specialized  term,  momentum  per  unit  charge,  that 

d.-sciihes  the  curvature  of  the  path  of  a  charged  particle  traversing  the  geo¬ 
magnetic  field.  It  is  a  unit  independent  of  particle  species  or  nuclear  com¬ 
position.  which  caii  he  converted  to  an  energy  for  any  specified  charged  nuclei. 
The  cutoff  rigidity  of  any  geographic  location  is.  a  function  of  both  the  zenith 
and  azimuth  angles  of  arrival.  At  the  surface  of  the  earth  and  in  the  vertic¬ 
al  direction,  the  cutoff  rigidity  has  a  value  of  ~13  to  -'18  GV  at  the  cosmic 
r  r.  y  equator  and  is  theoretically  zero  at  the  magnetic  poles.  In  actual  prac¬ 
tice  this  a* < magnetic  cutoff  is,  in  general,  not  sharp  and  in  more  technical 
leriis  we  i 'upper  computed  cutoff'  a  rigidity  above  which  all  particles  are 
allowed  in  the  s;  -.-cif  ied  direction,  'lowest  computed  cutoff',  a  rigidity  he  low 
which  no  par  Lie!**  has  tret,  computed  to  arrive  from  the  specified  direction,  and 
'•n  'effective  < -ut  of  f '  accounting  for  the  relative  transparency  of  the  cosmic 
c  a/  p  rr."  hr.:  (ih-.-:-  c -t  a.'. ,  IrtdO  . 


ih~  go’  ir‘r.i !  euu'l  ion  of  paiticle  motion  in  the  magnetic  field  does  not 
•solu-  ion  in  cUvd  form  even  in  a  simple*  dipole  field.  To  determ.’  r.  ■ 
igidit iok  are  a1  lev'-d  at  a  specified  geographical  location  and  direction 
cle  arrival,  it  :s.  newossa*.  y  to  pm  £«..:*■•!  ih-faj  1*  J  and  extensive  numeric- 
uhg  ions  of  cosmic  r-».  y  tiajc-cterio.-.  in  a  mat  hectical  model  of  t  h  • 
:,v 't;c  f  i-dd.  To  accural*1] v  detei  mjn.-.  the  cutoff  rigidity  of  a 
•e.  ion  on  i!.e  >■*. ict h  in  o  sv-cifu-d  direction,  cosmic  may  tra>*eU>:  - 

V*  J  ....  •  •  .  :  ,  n  aid  '  i . • 


•  «  *.  •_  *«  —  *■  _  »  ^  - 1  «  f  *  /LA-..  * 
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Estimates  of  cutoff  rigidities  obtained  from  trajectory  calculations 
vere  first  published  by  Freon  and  McCracken  (1962).  A  considerable  number  of 
such  cutoff  riaidities  have  subsequently  been  calculated  (Shea  and  Smart.  1962. 
and  references  therein). 

The  aim  of  these  calculations  has  generally  been  to  obtain  a  parameter 
suitable  for  the  ordering  of  data  from  cosmic  ray  detectors,  either  mobile  or 
stationary,  on  or  close  to  the  surface  of  the  earth.  The  cutoffs  obtained 
should  refer  to  the  epoch  in  which  the  data  were  obtained,  and  for  this  reason, 
and  also  for  the  purely  practical  reason  of  availability  of  field  models,  cut¬ 
offs  have  been  calculated  at  one  time  or  another  with  most  of  the  field  models 
which  have  been  available,  including  those  of  Finch  and  Leaton  (1957).  Jensen 
and  Cam  (1962).  and  1GRF  models  for  1965.0  (IAGA  Commission  2.  Working  Group 
4.  1969).  1970.0.  1975,0.  and  1975.0  projected  forward  to  1980.0  (IAGA  Division 
1.  Study  Group.  1976)  and  1980.0  (Peddie.  1982). 

2j_ _ METHOD 

To  specify  the  exact  path  of  a  charged  particle  in  the  geomagnetic 
field  it  is  necessary  to  numerically  integrate  the  differential  equation  of 
motion 


(r  =  5  t  *  b) 


a'ono  the  entire  trajectory  from  its  commencement  point.  In  principle  any 
ical  integration  technique  may  be  used.  Due  to  the  complex  trmectories 


(  I  ■«  Cl*  \  (l,:- 


. ic*  ray  p articles  with  which  we  are  concerned .  we  have 


si  andarJ.’  on  the  we  of  the  fourth  order  Punge-Kutta  integration  terhniuue 


2  •  .  ,  r  n  *  L -  .  •  *  »  -  >  *  •  *  *  J -  •  *  -  -  J  •»  -  «  _  A  »  .  n.kJk.lA  ,  - 


.->3  originally  deso  ih  .d  by  McC rack-:n  et  «?/. ,  1%2. 


For  prurti  c<<  1  reasons 


f  »,  * 


calculate  the  tra  ieetoi  y  of  a  negatively  char. god  particle  moving  outward  from 
the  earth  <<t  n  port  u.-uImv  direction  from  a  rp-.-c i f  j c  lor.. it. ion.  Such  «"i  trajec' 

cry  ic  identical  to  that  of  a  positively  charged  particle  of  equal  rigidity 
approociung  the  location  from  the  so  mo  direction,  although  traversed  in  the 
opposi  to  u  i  root  i  on . 

Each  of  the  cutoff  ricjidity  values  was.  calculated  in  the  same  mariner  by 
initiating  the  cosmic  ray  trajectory  at  the  top  of  the  sensible  atmosphere  (for 
this  purpose  an  altitude  of  20  km  above  the  surface  of  the  earth)  in  the 
vertical  direction.  The  calculations  were  continued  until  either  access  to 
the  interplanetary  medium  was  assured  (the  trajectory  extended  to  a  distance  of 
more  than  25  earth  radii.)  or-  the  trajectory  was  found  to  be  forbidden.  For¬ 
bidden  trajectories  were  divided  into  two  groups,  those  which  intersected  the 
solid  earth  (called  re-entrant  trajectories)  arid  those  for  which  no  solution 
could  be  obtained  within  a  reasonable  number  of  iterations  arbitrarily  selected 
at  200,000.  Trajectories  of  the*  latter  typo  are  generally  found  only  at  loca¬ 
tions  at  high  geomagnetic  latitudes.  They  do  not  form  a  significant  fraction 
of  the  results  reported  here. 


il0BLP...rciP...f>r‘  VERTICAL  CUTOFF  RlGIblTY  V?  1,1155 


.active  vertical  cutoff  rigidities  calculated  for  a  world  grit 


(fdi-'.i  and  -art .  1 9  -  J  util  icing  Uv*  lO^1  order  and  degree  I  GRP  mode]  for  bpoi  1; 


15r;0.0  (reddio,  19AZ)  are  illustrated  in  Figure  1.  Actual  cutoffs  were 
calculate  i  lor  an  altitude  of  20  km  for  ini. reiu-n'.s  5*  in  latitude  and  15*  in 
*  ■>  • •  0  . ,  b :!  ; I  i .  t  •  . .  Pu!  I  i-v.u  i.d.e  i nf  *  i  ; -ul.nl  i on  wu. 

ns-.-.l  to  d-d  rune  cstotfs  on  a  1*  !..y  1*  grid  fi . which  the  plcdUd  .,v..5 


wore  obtained.  Significant  differences  have  been  found  between  the  effective 
vertical  cutoff  rigidities  calculated  for  Epoch  1965.0  and  those  calculated  for 
Epoch  1980.0  Changes  of  the  order  of  ±0.10  GV  were  calculated  for  many  areas, 
of  the  world.  Major-  increases  were  calculated  for  the  North  Atlantic  Ocean 
area,  aTid  major  decreases  for  the  South  Atlantic  Ocean  and  South  American 
areas,  as  illustrated  in  Figure  2. 


4_, _ THE  COSMIC  RAY  HOP  ATOP 

When  it  became  practical  to  calculate  vertical  cutoff  rigidities  by  the 
trajectory-tracing  technique  it  also  became  feasible  to  locate  the  maximum  in 
the  vertical  cutoff  rigidity  as  a  function  of  latitude  at  a  given  longitude, 
thereby  determining  the  location  of  the  cosmic  ray  equator  by  theoretical  meth¬ 
ods  (Shea. 1969).  We  have  determined  the  cosmic  ray  equator  for  Epoch  1980.0 
in  the  following  manner.  Vertical  cutoff  rigidities  were  calculated  at  inter¬ 
vals  1°  in  latitude  and  5*  in  longitude  in  the  region  of  the  cosmic  lav  equat¬ 
or.  These  calculations  were  made  for  approximately  10  discrete  latitudes  at 
each  longitude  centred  around  the  location  of  the  equator  determined  for  pre¬ 
vious  Epochs  (Shea  and  Smart. 1975).  The  maximum  vertical  cutoff  rigidity  for 
each  longitude  was  determined  by  a  least  squares  fit  to  the  vertical  cutoff. 


i j(ii  lilies  along  each  longitudinal  meridian,  with  the  loci  of  th  t-j  5-*,o  points  d  o — 
fit.-:/:  the  cosmic  rav  equator.  The  location  of  the  cosmic  ray  equator  for 
lf'-;0.0  and  that  determined  in  a  similar  manner  for  epoch  1955.0  arc  shorn 
m  i'lotue  3.  The  cosmic  ray  equator  for  Epoch  1980.0  has  shifted  no.  /hoard  by 
r  ’■  ’tm  1*  in  latitude  between  longitudes  280 'E  to  360*  E  with  a  r.va;  shift 
t-r  'lwih'ly  i,-.-  t  e  than  4’  in  latitude  between  longitude  X  31G*F.  nr.J  3?5T.. 


in  the  cosmic  ray  equator  m  the 


*  Forth  American  arid  Atlantic  Ocean  areas. 


s'  i  rut  t  the 


This  shift  in  the  location  of  the  cosmic  ray  equator  is  consistent  with 
the  results  of  ffporrc*  and  Pomerantz  (1970)  who  reported  finding  a  similar  shift 
bv  experimental  measurements  along  longitude  346*.  Other  experimental  measure¬ 
ments  conducted  over  the  pest  45  years,  mostly  in  the  Pacific  Ocean  area,  are 
consistent  with  the  results  shown  in  Figure  3.  that  there  has  been  no  signific¬ 
ant  variation  in  the  position  of  the  cosmic  ray  equator  in  those  longitudes. 

5^ _ TPA3ECT0PV'  CALCULATION?  FOP  SATELLITE  ALTITUDES 

In  addition  to  the  calculations  reported  in  previous  sections  we  have 
also  computed  a  very  extensive  set  of  cutoff  rigidities  for  points  on  the 
shell  traversed  by  an  earth  satellite  in  a  circular  geocentric  orbit  at  a  mean 
altitude  of  400  km.  A  total  of  67  directions  of  particle  arrival  were 
considered  at  each  of  132  locations  distributed  each  10*  in  latitude  and  30*  in 
longitude  around  the  orbital  surface,  extending  in  latitude  between  50°N  and 
50’S.  Some  initial  results  have  been  reported  (Humble  et  ^<(..1979)  and  a  more 
detailed  report  is  currently  being  prepared.  The  calculations  were  commenced 
in  1977.  initially  in  connection  with  experiments  to  be  flown  on  board  the 
HEA0--C  spacecraft.  A.  field  model  for  epoch  1980.0  was  required,  and  the  IGPF 
1975.0  model  was  chosen,  extrapolated  forward  to  1980.0  by  use  of  the 
associated  secular  drift  coefficients  (LAGA.  1976.  co.cit.).  This  field  model 
is  of  8'^'  order  and  degree. 

Following  the  adoption  of  the  present  IGRF  for  1980.0  (Peddie . 1982)  we 
carried  cut  some  sample  tests  to  determine,  in  a  qualitative  fashion,  the 
effect  on  our  results  of  changing  from  the  model  previously  used  to  the  adopted 
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1980.0  model.  Due  to  the  amount  of  computer  time  involved  ue  were  anxious  to 
avoid  a  more  comorehensive  comparison  unless  such  were  forced  upon  us  by  the 
results  of  the  pilot  study. 


"ft  further  complication  exists,  induced  by  the  decision  to  define  tho 
adopted  1980.0  model  at  10^  order  and  degree.  Trajectory  integration  uses 
large  amounts  of  computer  time,  a  considerable  fraction  of  which  is  utilised  in 
evaluation  of  the  vector  components  of  the  geomagnetic  field.  The  integration 
technique  used  (section  ?.)  requires  these  magnetic  field  vectors  to  be  computed 
four  times  per  integration  step.  The  computer  time  required  for  each  computa¬ 
tion  of  vectors  increases  by  approximately  the  square  of  the  degree  and  order 
of  the  field  model  used,  and  the  change  from  8th  to  10th  order  and  degree 
increases  the  computation  time  for  a  given  trajectory  by  about  25%.  (That 
this  is  less  than  the  expected  value  is  partly  accounted  for  by  our  dropping 
the  higher  order  terms  in  the  expansion  when  they  become  insignificant  at 
increasing  altitude.)  We  have  therefore  also  undertaken  a  qualitative  survey 
on  the  effect  of  truncating  the  adopted  model  at  8th  order  and  degree. 


In  this  study  we  have  compared  cosmic  ray  cutoff  rigidities  calculated 
for  66  locations  on  our  world-wide  gr id  using 


(1)  the  8  order  and  degree  IGRF  1975.0  extrapolated  forward  to  1980.0 


(2) 


(abbreviated  as  X08  in  the  following  discussion). 

the  10  order  and  degree  IGRF  1980.0  adopted  field  model  truncated  to 
8th  order  and  degree  (abbreviated  as  AOS),  and 

r°  Ul(i  a-toptf-el  10th  order  and  degree*  IOT  1880.0  (Feddie .  1982)  abbreviated 
as  AID, 

Thfc  results  are  tabulated  m  table  1. 


■*«  1  ilhAi 


*  '  ^  »  V  •  »  •  «  ’  k  *  k  •  k  * 

» .  .... 


Inspect. ion  of  table  1  shows  significant  diff ercires  between  the  result 


obtained  uv.  ini  t  ho  various  models.  The  rigidity  intervals  used  were  Vi 
throuni K’ul .  abnormally  coarse  compared  with  many  invest j cud-ions,  but  it  is 
clear  oven  at  Ibis  interval  that  the  choice  of  field  mode-),  is  critical  in  some 
circumstances. 


Another  set  of  calculations  display  sensitivity  to  choice  of  field 
model.  Charged  particles  are  able  to  reach  orbiting  satellites  from  general  1 
westerly  directions  at  zenith  angles  considerably  larger  than  that  of  the  loca 
earth  horizon.  Wc-  have  carried  out  calculations  to  determine  the  largest 
zenith  angles  with  which  such  particles  are  able  to  reach  a  satellite  orbiting 
at  400  kn;  (Humble  et  at, .1983).  A  sample  of  the  results  obtained  using  the 
adopted  IG"^1  order  field  is  shown  in  figure  4a.  In  figure  4b  we  show  the 
differences  from  the  A10  field  results  obtained  by  use  of  the  MO 3  and  AOS 
models.  The  majority  of  the  differences  are  2°,  equal  to  the  zenith  angle 
increment  used. 


The  implication  of  our  results  is  that  the  choice  of  field  model  is 
import. ant.  even  at  the  relatively  coarse  sampling  intervals  in  rigidity  and 
dJ lection  of  arrival  which  have  been  employed.  However,  it  should  be  noted 
tint  this  is  a  general  statement .  which  does  not  apply  to  all  trajectories. 
In  m:iny  cases  individual  trajectories  are  rather  similar  regardless  of  the 
field  modal  used.  Figure  shows  altitude  vs  longitude  sections  for  such  a 
trajectory  calculated,  for  the  same  initial  conditions,  in  all  three  field 


IlKr.t’l;-.. 


There  is  very  little  significant  difference  between  the  three 
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For  (  i't  >ii  i  ra  .••dories,  however,  the  choice  of  field  model  is 
ctu.  i  i).  Fie;: ■: .  .  f.  and  7  depict  this  situation  for  two  different  sets  of 
uni  nl  conditi*  .i  The  figures  show  totally  different  trajectories  in  each 
of  the  field  mod.-ds.  The  figures  are  presented  for  illustrative  purposes 
only,  ay  an  indication  of  whnt  differences  can  be  found,  and  it  sh.uld  not  be 
concluded  from  them  that  such  differences  would  necessarily  be  found  over  a 
range  of  other  initial  conditions.  We  note,  however,  that  these  trajectories 
pass  through  the  Atlantic  region  longitudes  to  which  reference  has  been  made  in 
section  3. 

6, _ DISCUSSION  AND  O'  iMCLUSIOH 

In  the  past  few  years  it  has  become  more  apparent  that  the  cosmic  ray 
trajectory  calculations  used  in  the  analyses  of  cosmic  ray  observations  must  be 
made  utilizing  geomagnetic  field  models  for  the  epoch  in  which  the  data  were 
obtained.  The  general  secular  trend  that  has  become  apparent  over  the  past  25 
years  is  sufficiently  large  to  be  detectable  in  precise  cosmic  ray  measure¬ 
ments.  Furthermore,  it  would  seem  that  the  IGRF  for  Epoch  1975,0  and  its 
associated  predicted  time  derivatives  is  particularly  ill-suited  for  the 
analysis  of  cosmic  ray  data  acquired  in  the  longitude  embracing  the  Zitluntic 
Ocean  region. 


As  shown 


in  Figure  2b,  major  changes  are  found  between  the  vertical 


on',  of  {  rig  *  dil  i»  z  calculated  for  the  world  grid  for  Enochs  1965.0  and  1 9  >'< .  fj  in 
the  ntJontic  Ocean  region.  Tins  is  the  same  longitudinal  region  where  we  have 
found  r::grw f i c  uii  difference.';,  between  cutoff  rigidities  computed  for  1980  using 
li;i  *  i  >  iv  5  •  :•  •  or- 1  !'"-i  mure  do*  i  nit  ive  1'vO  magnetic  field  mode  l:;. 

Ti:;:-.  .is  also  tb  •  •:>  ea  where  latitude  surveys  have  shown  the  greatest  deviation 
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from  th**  result:;  expected  oil  the  basis  of  model  calculations  (Mischko  ft  at., 
1977  or  1979;  Shea  et  at. ,  1 c'  8 1 )  and  where  there  is  a  major  latitudinal  shift  in 
the  cornu c  ray  eguAtor  over  a  25-year  period. 


-We  note  that  regions,  in  which  we  find  major  changes  with  time  in  cosmic 
ray  trajectories  arid  cutoff  rigidities  are  similar  to  those  regions  in  which 
the*  horizontal  component  of  the  magnetic  field  shows  temporal  variations 
(Peddie, 1982) .  Cosmic  ray  particles  respond  to  the  integral  of  V  *  B  over 
their  entire  path  length  in  the  geomagnetic  field.  Cosmic  ray  trajectories 
having  rigidities  near  the  cutoff  have  a  larger  fraction  of  their  trajectories 
at  lower  altitudes  and  are  therefore  more  likely  to  be  affected  by  the  higher 
order  components  that  contribute  to  the  total  geomagnetic  field  topology. 
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Table  1.  Vertical  upper  computed  cutoff  rigidities  calculated  using 
the  extrapolated  (X03)  and  truncated  adopted  (AO 8)  geomag¬ 
netic  field  models  expressed  as  a  percentage  difference  from 
the  cutoff  calculated  for  the  10^  order  and  degree  adopted 
(A10)  field. 
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1  Estimated  vertical  cutoff  rigities.  epoch  1980.0  The  contours 
are  in  units  of  GV. 

2  Contours  of  averaged  annual  change  of  estimated  vertical  cutoff 
rigidity  between  1985.0  and  1980.0 

3  Estimated  geographic  location  of  the  cosmic-rav  equator  for  epochs 
1955.0  and  1980.0. 

4a  Largest  2enith  angles  at  which  charged  primary  particles  are  able 
to  reach  a  satellite  orbiting  at  400  Jun  altitude.  IGRF  1980.0 
adapted  to  10^  order  field. 

4b  The  largest  accessible  zenith  angles  found  using  other  field 
models,  expressed  as  a  difference  from  the  results  shown  in 
Figure  4a.  The  abbreviations  used  are  defined  in  the  caption  to 
table  1. 

5  Trajectories  of  15.844  GV  particle  arriving  at  (lat.0*.  long. 120’. 
alt.400!:m)  from  zenith  128'.  azimuth  225*. 

6  Trajectories  of  9.828  GV  particle  arriving  at  (0*.  300*.  400km) 
from  (128*.  If  .  *),  Note  that  the  arrival  direction  is  accessible 
from  infinity  a  the  A03  and  X08  field  models  but  is  inaccessible 


(trajectory  forbidden)  in  the  A10  field. 


Tro jectorios  of  10.537  GV  particle  arriving  at  0*,  GO*,  400km) 
from  (102*,  225*).  The*  direction  of  arrival  is  forbidden  from 
infinity  in  both  the  AOS  and  A10  fields,  but  is  accessible  accord 
iny  to  the  earlier  (X03)  model. 
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Hits  Below  the  Belt  -  Cosmic  Ray  Direct  Access  to  Low 
Altitude  Satellites  from  Earthward  Directions 

J.  E.  HUMBLE 
Department  of  Physics 
University  of  Tasmania 
Hobart,  Tasmania  7001,  Australia 

D.  F.  SMART  and  M.  A.  SHEA 
Air  Force  Geophysics  Laboratory, 

Hanscom  AFB,  Bedford,  MA. ,  01731 

ABSTRACT 

We  have  previously  reported  finding  that  primary  cosmic  ray  particles 
are  able  to  reach  satellites  at  400  km  altitude  from  zenith  angles  up 
to  133°.  More  detailed  investigations  have  now  found  that  access  is 
possible  at  certain  locations  at  zenith  angles  up  to  150°,  only  30°  from 
the  nadir.  Access  from  zenith  angles  of  order  140°  to  144°  is  common  and 
occurs  at  all  the  latitudes  which  we  have  investigated  between  40°  North 
and  30°  South.  Particles  arrive  at  these  large  zenith  angles  from  azimuths 
on  the  equator  side  of  west  in  both  hemispheres. 

Presented  at  the  Fall  meeting  of  the  American  Geophysical  Union, 

San  Fransisco,  CA,  December  1982. 


The  Directions  from  which  Cosmic  Rays  may  Reach  Earth  Satellites 
by  HUMBLE.J.E.,  UNIVERSITY  OF  TASMANIA 

It  is  easy  to  feel  intuitively  that  primary  cosmic 
rays  are  able  to  reach  earth  satellites  from  all  direc¬ 
tions  above  the  horizon  and  will  be  unable  to  do  so  from 
directions  much  below  the  horizon.  The  falsity  of  the 
latter  assumption  is  discussed,  and  it  is  shown  that 
primary  particles  can  reach  satellites  at  400  km  altitude 
from  zenith  angles  up  to  130°.  At  an  altitude  of  1250 
km  primaries  may  gain  access  to  the  satellite  from  zenith 
angles  up  to  at  least  178°,  essentially  from  directly 
beneath  the  satellite. 

Presented  at  the  annual  meeting  of  the  Astronomical  Society  of  Australia, 
Sydney,  N.S.W.,  Australia,  May  1983. 


FIELD  MODEL  COMPARISONS:  COSMIC  RAY  DIFFERENCES  BETWEEN  THE  PREDICTED 

AND  ADOPTED  1980.0  FIELDS 


J.E. Humble 
Physics  Department 
University  of  Tasmania 
Hobart,  Tasmania,  Australia 


Cosmic  ray  trajectories  have  been  calculated  using  the  10th  order  definitive 
geomagnetic  reference  field  adopted  at  the  1981  1AGA  meeting  in  Edinburgh; 
the  same  field  arbitrarily  truncated  to  8th  order;  and  the  Interim  Geomagnetic 
Reference  Field  for  1980.0,  based  upon  the  published  1975.0  field  and  its 
associated  secular  drift  coefficients.  It  is  shown  that  in  the  majority  of 
cases  the  trajectories  obtained  using  each  field  model  are  rather  similar. 
However,  in  certain  circumstances  considerably  different  trajectories  result, 
particularly  if  the  location  concerned  is  near  to  the  South  Atlantic  Geomagentic 
anomoly.  The  use  of  the  different  field  models  does  not  significantly  alter 
the  largest  accessible  zenith  angles  calculated  for  some  sample  locations. 

Presented  at  the  XVIII  General  Assembly  of  the  International  Union  of  Geodesy 
and  Geophysics,  International  Association  for  Geomagnetism  and  Aeronomy 
meeting,  Hamburg,  Federal  Republic  of  Germany,  August  1983. 


COSMIC  RAY  ACCESS  TO  SATELLITES  FROM  LARGE  ZENITH  ANGLES 

J.E. Humble 

Department  of  Physics 
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and 

D.F. Smart  and  M. A. Shea 
Air  Force  Geophysics  Laboratory  (PHG) 
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ABSTRACT 

As  the  altitude  of  earth-orbiting  satellites  increases 
it  becomes  progressively  easier  for  primary  cosmic- 
ray  particles  to  gain  access  to  them  from  directions 
below  the  geometric  horizon.  Results  are  presented 
from  a  comprehensive  survey  of  the  largest  accessible 
zenith  angles  for  a  range  of  altitudes  and  geographic 
locations.  The  search  has  disclosed  that  primary 
particles  are  able  to  reach  a  satellite  at  an  altitude 
of  1250  km  from  zenith  angles  as  large  as  178°. 

Presented  at  the  18th  International  Cosmic  Ray  Conference,  Bangalore,  India, 
August  1983. 


Cosmic  Ray  Access  to  Spacecraft  from  Earthward  Directions: 
the  Role  of  the  Atmosphere 

J.E. Humble 

Physics  Department,  University  of  Tasmania 

Previous  calculations  on  this  topic  have  used  a  very  simple 
model  for  the  atmosphere.  The  model  assumed  the  atmos¬ 
phere  to  be  totally  transparent  to  cosmic  rays  at  altitudes 
above  30  km  and  totally  opaque  to  them  at  lower  altitudes. 

The  calculations  reported  here  use  a  realistic  model  atmos¬ 
phere.  It  is  shown  that  the  results  obtained  are  similar 
to  those  obtained  with  the  cruder  model. 

Presented  at  the  annual  meeting  of  the  Astronomical  Society  of  Australia, 
Coonabarabran,  N.S.U,  Australia,  May  1984. 
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